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A systematic dynamical system approach is applied to study the cosmology of anisotropic Bianchi 
I universes in which a vector field is assumed to operate on a disformal frame. This study yields a 
number of new fixed points, among which anisotropic scaling solutions. Within the simplifying as¬ 
sumption of (nearly) constant-slope potentials these are either not stable attractors, do not describe 
accelerating expansion or else they feature too large anisotropies to be compatible with observa¬ 
tions. Nonetheless, some solutions do have an appeal for cosmological applications in that isotropy 
is retained due to rapid oscillations of the vector field. 


I. INTRODUCTION 

The observed large-scale homogeneity and isotropy of the Universe are usually explained within the framework of 
cosmological inflation. This epoch of primordial exponential acceleration is traditionally assumed to be driven by a 
scalar field; one immediate consequence of this choice is that whatever primordial anisotropies might be present at 
that time, they are rapidly washed out as inflation proceeds, leaving the Universe with an order 10“ 5 initial scalar 
(and somewhat smaller tensor) seeds for the observed structures, germinating from the quantum fluctuactions of the 
inflaton scalar field itself (see [l] for a review). Scalar fields are also often invoked in models of late-time dark energy, 
according to which the present acceleration of the Universe is due to the presence of a tiny, smoothly distributed 
vacuum energy associated with the so called quintessence field (see, e.g., [2|). 

Despite their attractiveness and simplicity however, the existence of these cosmological scalar fields has not yet been 
established observationally, leaving the window open for other possibilities, among which the logical simplest option 
is vector fields. Vector fields in this context could either play the roles of cosmological prima donnas as either the 
actual inflaton or quintessence fields driving the early or the late acceleration, or as supporting acts to a scalar field, 
which may have non-trivial observable consequences. For instance, one motivation for entertaining such ideas is the 
apparent tension among several features in the observed cosmic microwave background radiation and the statistically 
isotropic standard model of the Universe, see e.g. Ref. HE] for recent studies. 

Of course, in any reasonable model of our Universe, isotropy (at least approximate) must be had in order to not 
clash with observations; this means that one has to consider non-minimal options for the vector fields configurations in 
one way or another. One could support isotropic expansion long enough with sufficiently finely tuned initial conditions 
as in the vector field inflation model of [5] , or invoke a regiment of vector fields in such a way that on average they 
do not manifest any preferred direction, as in the model put forward by EH3, see also EH2H- With non-minimal 
coupling to gravity or non-standard vector fields, such as in three-form inflation and dark energy models, a time-like 
(and thus isotropic) field would include non-trivial (yet observationally appropriate) cosmologies [IBfEijl . Also, one 
may arrange a “triad” of three vector fields in such a way that their vacuum configuration is isotropic m- Finally, 
rapidly oscillating vectors could naturally yield isotropic cosmology |28j . 

Here we investigate a novel class of scalar-vector cosmologies. In particular we are interested in cosmological 
dynamics that could yield the so-called “anisotropic hair”, i.e., vectors which survive the inflationary stretch. Such 
scenarios would manifest themselves as attractor solutions in which the isotropic scalar field dominates the energy 
budget of the Universe, but a small, constant amount of anisotropy is supported by subdominant vector fields. Such 
solutions have been found in non-minimally coupled scalar-vector models [29MT1 . see also [32H35] for studies of the 
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dynamics of anistropic vector cosmologies. Here we will study a system in which the non-minimal coupling is of the 
disformal type: the (generically massive) vector field is thus coupled to the gradients of the scalar field. We study a 
simple axisymmetric Bianchi I spacetime and its anisotropic cosmology, in which the direction of the space-like vector 
determines the preferred axis. 

The formalism we employ was recently introduced in Ref. l¥i'i . In that case we focussed on a particular extra- 
dimensional theory in which the vector matter lives on a moving brane |37l439j . giving rise to a disformal coupling 
of the vector to the scalar. Since in that simplest model we could not realise any “anisotropic hair” solution we 
are led to explore minimal tweaks to it with the hope of finding new more interesting cosmological solutions. The 
class of Lagrangians we study here are a “doubly-boosted” version of that studied in our previous work, thanks to 
an additional Lorenz boost factor (due to the brane movement). The tools for the phase space analysis developed in 
Ref. j36j . as anticipated in that work, are extremely handy for the present class of models, and this work is one more 
proof of their versatility. After briefly defining our double-boost class in the following section [TTJ we formulate the 
cosmological equations of motions in the language of dynamical systems using the formalism of Ref. |36| . The actual 
analysis is performed in section |TV} and we conclude in section [VT| 


II. THE ACTION 


The action for this analysis is given by 


S = 



R 

167tG 


+p(<t>,x) 


Sa- 


(1) 


The first piece includes standard Einstein gravity and the scalar with a lagrangian that depends on the field and its 
kinetic term X = —g a P(j) t0 ,<f) t p/2. For a canonic scalar field one has p = X — V, where V(cf>) is the scalar potential (p 
is in fact the “pressure” of the field in fluid description). In our notation a comma stands for partial derivatives. The 
second piece is a disformally-inspired Proca lagrangian for the vector field, which explicitly reads 


S A = d A x^g 


-2 


--F 2 

4 


-M 2 A 2 

2 


+ 


+ M 2 A a Ap ) 


( 2 ) 


where K = K{(j>) and M 2 = M 2 (</>). The field strength tensor is F a p = Ap, a — A a ,p as usual. The Lorentz-like factor 
is given by 

7 -2 = 1 — 2 KX . (3) 


This action can be obtained from a more compact form realisation of a system in which the vector field lives on a 
disformal metric of the kind 




(4) 


where C(cf >) and D(</>) are the conformal and disformal factors respectively. In this work in particular we study a 
variant of the simpler form which was investigated in !36| . In general one finds that I\ = D/C. 

The next step is to obtain explicit covariant field equations for this system, which are 

G, V =K 2 {T*+TP , (5) 

where the gravitational coupling constant is a 2 = 1 / 87 tG; the stress energy tensor for the scalar field 4> is 

T$ v = pg + p,x<t>,»<l>, v , ( 6 ) 

while that of the vector is 


T A _ p p 7 , 


7 -2 ( -\f 2 - ^M 2 A 2 ^j + (F 1CL F\ + M 2 A a A^ 


9nv 


K 


\</>A,»F 2 - <p’ a ^F m F„p + <\>^F m F v a - t^F^FJ - (t>, v <j> ,ol F ai F 


+ M 2 [7 2 A IJb A l , + K(j)^(j) tV A 2 + K {(j)^A v + (f> jV A^)j . 


(7) 


Now, since there is an explicit coupling term between the scalar field and the vector, the stress-energy tensors of the 
two are not separately conserved, but, as it must, + T'/ v ) = 0. 

We now set out to study these equations of motiong in a general case, where, in order to allow for the possibility 
that anisotropies can develop, we have to employ anisotropic spacetimes. 
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III. COSMOLOGY 

We choose to work with a Bianchi I axisymmetric metric, since it describes an anisotropic expansion with the 
anisotropy pointing in one particular direction (we choose the latter to be aligned with the z direction). The metric 
is 


ds 2 


-dt 2 + e 2a{t) 


e 2 ^ ( dx 2 + dy 2 ) 



( 8 ) 


Here a is the overall isotropic volume expansion rate (the Hubble parameter), and a is the anisotropic shear. We can 
project the covariant equations of motion for this metric with the result: 


0 

0 


0 

0 

0 


A + (d + 4 a) A + (l - A> 2 ) M 2 A , 

($ + 3 otifij p, x + 2 p,xx + P,<i>x4> 2 - P,<j> + 2KM 2 e~ 2a+ia AA<fi 


g-2a+4cr^2 


KM 2 (j> + (a + 4 &)^j + \k'M 2 ^ 2 + MM' (l + Kft 


3 cr + 9dd - e" 2 Q+ 4 CT i 2 + e~ 2a+4a A 2 M 2 (l - K<jA , 


6 d + 9d 2 + 9(i 2 + 3 p+ ^e~ 2a+4a A 2 - ^e~ 2oi+4a A 2 M 2 (l - K^ , 
3d 2 - 3d 2 -p,x ^ 2 +P~ \ e~ 2a+4,T A 2 - i e - 2a+4a A 2 M 2 (l + AA/ 2 ) . 


At this point we find it convenient to define the inverse of the Lorentz factor 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 


7 = 1/7 = y 1 — Kef ) 2 . 


(14) 


It is most suitable, for studying the dynamics of the system, to employ the e-folding number x = a as the time- 
evolution variable. Also, to be concrete, we choose the scalar field action as the DBI type 

After several manipulations we found that the best first-order variables which make the system the most transparent 
and easy to follow are 


1 


x = 


^37(7+1) 


, y = 


k 2 V 

3d^’ 


u = e~ a+2a 


yj2 - 7 2 M . e 

v^d 


— a+2cr 

kA , v = -— kA 


V6 


M 

,x, M = — . (15) 

a 


Notice that A4 simply tells us how M evolves compared to the overall expansion rate, and does not encode any direct 
dynamic effect upon the background geometry. The derivatives of the metric are described by the shear and the 
slow-roll parameter 


E = a x , e = — d/d 2 


The Friedmann equation (13) then beautifully reads 

1 - E 2 = x 2 


I Z 1 Z 1 z 

+ y +u + v 


and the second Friedmann equation (12) constrains e, almost as beautifully, as 

2e = 3(E 2 + 1) + 37 a ; 2 - 3 y 2 - (2q 2 - l)” 1 ^ 2 + ' 


(16) 


(17) 


(18) 


To close the system we choose to work with the simplest slopes for the three (/-dependent functions, that is, we assume 
them to be approximately constant, which is equivalent to assuming a nearly exponential form for these quantities. 
Therefore, we can define the, approximately time-independent, slope parameters 


Ay = 



A k = \ hz 


3 K. 


2 kK 


Am 


= 2 


V 2 kM ' 


(19) 
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All ingredients are in place for the dynamical equations, which can be cast as 


where 


) 2 X, M = (1 + 27)(1 - 7 ) [v^AjfX - e] x + T, 

( 20 ) 

y,u = (V*X v x + e)y, 

( 21 ) 

u, m = (^yfk\MX + e + 2£ — 1 + '"^ 2 ^ u + \/^ ~ l 2 Mv , 

( 22 ) 

7 2 

v jv = (e — 2E — 2) v -. Mu , 

\/2-7 2 

(23) 

— (e 3) E + 2v + ~ u , 

7 z — z 

(24) 

hA. H - , 

(25) 

77, AA = ~ {(7+1) [V*X K x 2 ex] + T} ; 

(26) 


notice that we have eliminated I\ thanks to 3 x 2 K = 7 — 1 = (1 — 7 )/ 7 - We have kept aside the second derivative of 
the scalar field, which only appears in the combination 


T = 


_ / 7 + 1 

37 


ft'/’.A/'Af 


(27) 


{\/*z 2 \l A ( 2 A K x 2 - 2 A v y 2 - 2 (A k + A M )u 2 ) 


2(7 — 1)7 2 m 2 + ( 7 2 — 2)x 2 
7^A kx 2 -f- 7 2 (— 5Xkx 2 -f~ 4Ayy 2 A 2(A/v + 2 Am)'^ 2 ) + 2^yXjex~ + 2A/e ; u] 

— (7 + \)x ( 37 2 — e)( 7 2 — 2)x 2 + 47 2 (7 — 1)\/ 2 — 7 2 M.uv + 27 2 (7 — 1) (4£ + 1 — e) m 2 J | ; 


this combination is manifestly finite in the two extreme limits of 7 —> 1 and 7 —> 0 , which are those we will focus on 
in the subsequent analysis. Finally, we use the shorthand yfk = ^ 7(7 + l)/2. 

We can use the Friedmann constraint to solve away one of the variables in the system, which we choose to 
be y assuming it is always positive 1 ; the phase space is therefore a four dimensional one. As anticipated above, we 
study the system in the so-called ultra-relativistic (7 = 0 ) and non-relativistic (7 = 1 ) regimes. 

It is interesting to notice a peculiar feature of this system: the coupling of the scalar and the vector is turned off 
when the mass vanishes M = 0, that is, when the vector field is conformally coupled to g it also decouples altogether 
from the scalar held (j) in the physical metric g Nonetheless, this does not mean that the dynamics of the system 
are completely decoupled if M. —>■ 0 , since M quantifies the “weight” of the vector held with respect to the expansion 
rate, but there is additional dynamical information about the vector potential in u 2 : if the vector held is very light 
in units of the expansion it could anyhow affect the expansion and the shear. 


IV. CLASSIFYING THE FIXED POINTS 
A. Ultra-relativistic regime 


In the 7 —> 0 case we obtain a very simple closed system: 

(x 2 ) jAf = (2e-3)x 2 - (l + 4E)zt 2 - 2y/2Muv , (28) 

it iAf = (e + 2E — l)rt + , (29) 

v = (e — 2E — 2)v , (30) 

E iAf = (e — 3)E + 2v 2 , (31) 

-Mw = eM. (32) 


1 


Negative energy density potentials V < 0 could be allowed by extending the analysis to allow imaginary y. 
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Quite remarkably, there is no more reference to the three slopes Ay, A k 1 and A m- this means that the structure of the 
scalar functions C(<f> ), D{(f>) 7 and M{(f> ) do not influence the dynamics directly in this limit; nonetheless, the stability 
properties do depend on Xx since it is this slope which tells us whether 7 —> 0 or not. The fixed points are in this 
case identical to the warped brane scenario: 1) de Sitter expansion with y = ±1 for any XA driven by the potential 
of the scalar field; 2) kinetic domination with x = ±1 (and everything else set to zero); 3) anisotropic solutions with 
£ = ±1. Thus, either there is no anisotropy, or the anisotropy is very large. 

B. Nonrelativistic regime 

Let us move on to the limit where 7 —► 1. It is clear at first glance that this case is much more interesting: 


£,jv = (e - 3)x - Am « 2 - A vy 2 , (33) 

u = (e T 2£ 1 T A mx)u T XAv , (34) 

u jjV - = (e — 2£ — 2)v — A 4u , (35) 

S )Af = (e — 3)£ — 2u 2 + 2v 2 , (36) 

XA^x — (e f- \mx)M . (37) 


In order to find the fixed points of the system we begin with M. —> 0, i.e., we look at solutions for which the vector 
is light with respect to the expansion rate; this trivially solves the last equation. We were also able to obtain some 
fixed points with arbitrary A 4 but only for very special values of the coefficients (A^, Am) which we list together with 
their XA = 0 relatives. 

In studying the stability of the fixed points we follow the definitions of Lyapunov as in our earlier work on this 
subjet, see [36] : 

Once again, in order to check whether 7 = 1 is a stable solution we need to expand around this value: the coefficient 

of the linear term in (1 — 7 ) > 0 , whose sign we need to determine, is given by: 

(! - 7 ),w = -2 (3 - \ K x + A v y 2 /x + X m u 2 /x) (1 - 7 ). (38) 

There are five families of fixed points, which are described in turn below. For each fixed point, we first specify the 

solution, then write down the conditions for its existence and stability, and then briefly comment upon its potential 
applications in cosmology. 


1. Isotropic scalar field domination 


fz = -Ay/3, y=\J 1 - Ay/9 , 

| a = 0 , v = 0 , 

\ E = 0 , and e = Ay/3 . 


( 39 ) 
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Ad = 0 

Existence 

Ay < 3 

Stable 

0 < Ay < y/6 and A m > Ay 
— y/6 < Ay < 0 and Am < Ay 

Senristable 

Am = Ay , Ay 1=3, Ay|=\/6, Am = Ay —3/Ay 


0 

Existence 

A m — Ay or Ay = 0 

Senristable 

6 Ad < Ay - 3 

Lyapunov 

6 Ad > Ay - 3 


7 —> 1 

A a' > 0 and either Ay < 0 or Ay > — Ar- 
A a < 0 and either Ay > 0 or Ay < —Xk 

> Comments - This is a valid accelerating solution. Not the fixed point itself but its stability properties are affected 
by the vector field. 


2. Kinetic driven anisotropic universe 


! x = x, y = 0, 

u = 0, v = 0, (40) 

E = sjl — x 2 , and e = 3 . 

> Comments - This point with Ad = 0 is not stable: there is always one eigenvalue which has positive real part. 
Also the Ad 0 fixed point defined by x = —3/A m and \Xm\ > 3 are not stable. 


3. Anisotropic scaling solution 1 


6Ay 

X = ~X 2 v + 12 
u = 0, 

v _ 2X 2 V - 12 

Ay + 12 


3x/24- 2X 2 
V= X 2 v + 12 ’ 

, \J (Ay — 12)(18 — 3Ay) 
V X 2 y + 12 


< 


(41) 
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A4 = 0 

Existence 

y/6 < | Ay | < 2x/3 

Stable 

a/ 6 < Ay < 2 a/3 and Am > Ay 
—2 a/ 3 < Ay < — a/6 and Am < Ay 

Semistable 

Ay = 2 a/ 3, Ay = a/6 

a/ 6 < Ay < 6 a/3/17 and either Am = Ay or Am = 3(Ay/2 — 2/Ay) 

Lyapunov 

6a/ 3/17 < Ay < 2 a/3 and either Am = Ay or Am = 3(Ay/2 — 2/Ay) 


A4 7^0 

Existence 

Am = Ay and either Ay = 2a/3 or Ay = a /6 

Semistable 

Ay = a /6 and A4 < 1/2 

Lyapunov 

Ay = a /6 and A4 > 1/2 


7 —7 1 

A k > 0 and either Ay < 0 or Ay > —A k 
A k < 0 and either Ay > 0 or Ay < —A k 

> Comments - This solution is stable in some parameter range, however it cannot combine accelerating expansion 
with small shear since one always has £ = —l + e/2. 

4- Anisotropic scaling solution 2 


A m ~ 2Ay 


1 i/9Ay — GAyfAy — 36 

2 A m ~ 2Ay 
v _ 1 Ay — 2A M 

2 A M ~ 2Ay 


1 \/6A| f — 9 Ay/Ay + 36 
^ 2 A M — 2Ay 

_ 1 a/ 6Ay — 6A^- + 3 Am Ay — 36 

2 Am — 2Ay 

3Ay 

and e = — r~ • 

2Ay — Am 


(42) 


Ad = 0 

Existence 

(i) |Ay| = 4 a/ 6/17 and |A M | = V 6 / 17 

(ii) Ay = 2a/ 3 and either Am = 2a/ 3 or Am < a/ 3 
(hi) (3Ay + a/9A^ - 96) < 4A M < (Ay + a/17A^ - 96) 

and either Ay < —2a/ 3 or 4 a/2/3 < Ay < 2a/3 

(iv) (Ay - a/17A^ - 96) < 4A m < (Ay + y/UX^ - 96) 

and 4 a/ 6/17 < |Ay| < 4^/2/3 

(v) (Ay - y/17Xy - 96) < 4A m < (3Ay - a/9A \ - 96) 

and either Ay > 2a/ 3 or 4 a/2/3 < Ay < 2a/3 

Semistable 

Am = Ay with | Ay | = 2a/3 


Ad 7^0 

Existence 

Am = Ay and | Ay | = 2a/3 


> Comments - This hxed point is unstable everywhere but for the very specific choice Ad = 0, Am = Ay with 
|Ay | = 2\/3, for which u = v = 0 and £ = 1/2, see above. 





5. Anisotropic scaling solution 3 


9Ay 


u = ± 


S = 


2A \ + 12 ’ 
3^9 ~ 3Ay 
2Ay + 12 
4Ay — 3 


y = o, 


V = =b 


^54 + 18Ay - 12Ay 


v 


2Ay+12’ 


and e = 


2Ay + 12 
6 A y + 9 
A v + 6 


M = 0 

Existence 

Am = 2Ay and | Ay | < \/ 3, or | Ay | = -\/3 


M ^0 

Existence 

| Ay | = y/3 and A m = 2Ay 


E> Comments - Nowhere in its domain is this solution stable. 


6. Anisotropic scaling solution 4 


(43) 


__ 3(Am — 2Ay) _ 

3A ^ — 5Am Ay + 2Ay + 12 
; \Z(3A^ - 4A M Ay + 12)(3A^ - 3A M Ay ^9) 

3A^ — 5A M Ay + 2Ay + 12 
_ 2(Ay — Am Ay — 3) 

, 3A^ — 5 A M Ay + 2Ay + 12 


\/(3A m — 4A M Ay + 12)(3A^ — 3XmXv + 18) 
3A|, / -5A M Ay + 2A|. + 12 


v = 0 , 


and (i) e = 3 , or (ii) e 


9Ay 

10Ay + 12 ' 


M = 0 

Existence 

(i) Ay = 3(Am/4 + 1/Am) 

(ii) Ay = -Am and |A M | > \f%j% 



Existence 

Am = Ay and |Am| = 2-\/3 


> Comments ■ Once again, there is no stable domain. 


(44) 


V. NUMERICAL EXAMPLES 

Let us gain some intuition about some of these solutions by looking at the numerical evolution of the system. The 
initial conditions are Xq = yo = 1/4 and Uq = Vq = 1/2. 

We begin with one of our attractor solutions, the Lyapunov stable point Ay = Xk = Am = —3, which we follow for 
7 o = 1/3 and A4o = 1/10, see Figure]!] Here 7 quickly attains its non-relativistic attractor, while the vector and the 
scalar quantities all oscillate with decaying amplitudes: the limit for t —> 00 does not exist (due to the oscillations) 
but the point is nevertheless a Lyapunov attractor, as all quantities remain bounded, including the non-zero shear S. 

In the following three examples we set Am = — A_r/ 2 2 . In Figure [ 2 ] we plot the model Ay = —5, A k = 8 . We 
start from 70 = 99/100 and a heavy Ado = 10. In this case, the asymptotic future attractor will instead be an 


2 This choice is reminiscent of a stringy interpretation of this model in terms of braneworlds. If h = D/C ~ exp (2 kA</>)> then C ~ 1/D 
exp (—kA</>). 
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FIG. 1: The model Xv = Xk = A m = —3, with 70 = 1/3 and Mo = 1/10. Left panel: The scalar field evolution. The dashed 
red line is x, the solid orange line is y. The 7 is the thick dotted blue line and the thick purple line is the total equation of 
state. Middle panel: The vector field evolution. The dashed red line is v, the solid orange line is v. The shear E is the thick 
purple line. Right panel: The phase space evolution in the ( u , v ) plane. 





FIG. 2: The model Av = —5, A k = 8 , A m = —4, with 70 = 99/100 and Mo = 10. Left panel: The scalar field evolution. 
The dashed red line is x, the solid orange line is y. The 7 is the thick dotted blue line and the thick purple line is the total 
equation of state. Middle panel: The vector field evolution. The dashed red line is v, the solid orange line is v. The shear E 
is the thick purple line. Right panel: The phase space evolution in the ( u , v) plane. 


ultra-relativistic (7 —> 0) with contributions to the expansion rate from both the vector and the scalar held. We see 
that the vector held and the scalar held oscillate before settling into the attractor: this is an effect of the large M, 
which will grow indefinitely. Strictly speaking this is not a pure attractor because of the unbounded behaviour of M , 
but as we already mentioned, there is no dynamics attached to this variable. 

In Figure[3]we show the case of Ay = 5, A k = 3. Here the asymptotic future attractor will again be a non-relativistic 
(7 —> 1) solution. The vector will be strongly oscillating, so that neither u nor v are constants. However, their total 



FIG. 3: The model Xv = 5, A k = 3, Am = —3/2, with 70 = 99/100 and Mo = 10“'. Left panel: The scalar field evolution. 
The dashed red line is x, the solid orange line is y. The 7 is the thick dotted blue line and the thick purple line is the total 
equation of state. Middle panel: The vector field evolution. The dashed red line is v, the solid orange line is v. The shear E 
is the thick purple line. Right panel: The phase space evolution in the (u , v) plane. 
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V 





FIG. 4: The model Ay = —A k = 2A m = —4, with 70 = 99/100 and A4o = 1/30. Left panel: The scalar field evolution. The 
dashed red line is x, the solid orange line is y. The 7 is the thick dotted blue line and the thick purple line is the total equation 
of state. Middle panel: The vector field evolution. The dashed red line is v, the solid orange line is v. The shear S is the 
thick purple line. Right panel: The phase space evolution in the ( u , v ) plane. 


contribution to the expansion rate, u 2 + v 2 , will be approximately a constant, and will scale with the scalar field 
contributions x 2 and y 2 with on average, a nearly dust-like equation of state. This is an interesting solution as it 
shows that a spatial vector field is compatible with the FLRW symmetries, and can mimic a matter-like component. 
Note that the shear decays in the evolution even though the vector field contribution does not: this is line with the 
isotropy theorem for cosmic vector fields put forward by in [28| . 

Finally, in Figure [4] we plot the special case Ay = —A k = 4, with again 70 = 99/100 but now non-negligible 
A4o = 1/30. When the A k and Ay have equal magnitude but the opposite sign, fixed points with finite 7 appear for 
the scalar field when it is dominated by its kinetic term. Similar solutions for the DBI scalar were reported in the 
absence of the vector m ■ Once again notice how the overall shear decays, FLRW is recovered, albeit the contributions 
from the vector oscillate strongly and are not tamed. The equation of state also strongly oscillates, and its final value 
depends ultimately on the choice of Ay. 


VI. CONCLUSIONS 

In this work we studied the cosmological dynamics of a particular vector-tensor model of gravity which is a simple 
modification and extension of the simplest disformal vector-tensor theory, which arises in presense of a double boost 
factor for the vector sector. Because of the presence of a single vector field, which can develop a non-trivial vacuum 
expectation value, we worked with an anisotropic Bianchi I spacetime, which includes anisotropic shear. We were 
able to reduce the system to a first-order one which can be studied with usual dynamical systems techniques. The 
equations of motion follow the evolution of a massive Abelian vector field, a generic real scalar field, and gravity, all 
of which were parametrised in such a way that the relatively cumbersome system becomes fairly manageable. We did 
obtain a number of new fixed points, whose existence and stability properties we then studied systematically. For the 
sake of concreteness, we picked exponential forms for both the scalar and vector potential and their coupling, that is, 
in first approximation, we kept the slopes Ay, A m, and A k as constant. 

The outcome of our analysis does share several aspects with the single-boost warped brane case we investigated 
in our previous work, in the sense that, also with the modification presented in this work, the shear, which could in 
principle be supported by anisotropic stress of the spatial vector field, decays and the metric isotropises towards flat 
FLRW. Compared to previous work however we were able to obtain even more anisotropic fixed points, including 
scaling solutions; they however still seem not to be of immediate cosmological relevance because either they are not 
attractors, or it is not possible to have simultaneously a slowly rolling vacuum energy driving the expansion, and 
small shear. Curiously, one can show that for all the fixed points we obtained, there exists an approximate relation 
e ~ E + 0(1): this means that either e < 1, or E C 1, but not both. 

This however does not mean that there is no cosmological application for this theory. Indeed, as illustrated in our 
numerical examples, even in the absence of shear, the single vector can be compatible with the FLRW symmetries due 
to its very rapid oscillations, as predicated in [25] , Once the vector oscillated very rapidly its equation of state can 
be, on average, like that of dust, which in turn implies that such models can provide a disformally coupled alternative 
to dark matter. 
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